Search for: All records

Abstract Pulsar halos are regions around middle-aged pulsars extending out to tens of parsecs. The large extent of the halos and well-defined central cosmic-ray accelerators make this new class of Galactic sources an ideal laboratory for studying cosmic-ray transport. LHAASO J0621+3755 is a candidate pulsar halo associated with the middle-aged gamma-ray pulsar PSR J0622+3749. We observed LHAASO J0621+3755 with VERITAS and XMM-Newton in the TeV and X-ray bands, respectively. For this work, we developed a novel background estimation technique for imaging atmospheric Cherenkov telescope observations of such extended sources. No halo emission was detected with VERITAS (0.3–10 TeV) or XMM-Newton (2–7 keV) within 1^∘and $10^{\prime }$around PSR J0622+3749, respectively. Combined with the LHAASO Kilometer Square Array (KM2A) and Fermi-LAT data, VERITAS flux upper limits establish a spectral break at  ∼1–10 TeV, a unique feature compared with Geminga, the most studied pulsar halo. We model the gamma-ray spectrum and LHAASO-KM2A surface brightness as inverse Compton emission and find suppressed diffusion around the pulsar, similar to Geminga. A smaller diffusion suppression zone and harder electron injection spectrum than Geminga are necessary to reproduce the spectral cutoff. A magnetic field ≤1μG is required by our XMM-Newton observation and synchrotron spectral modeling, consistent with Geminga. Our findings support slower diffusion and lower magnetic field around pulsar halos than the Galactic averages, hinting at magnetohydrodynamic turbulence around pulsars. Additionally, we report the detection of an X-ray point source spatially coincident with PSR J0622+3749, whose periodicity is consistent with the gamma-ray spin period of 333.2 ms. The soft spectrum of this source suggests a thermal origin. 

Abstract Assuming Galactic cosmic rays originate in supernovae and the winds of massive stars, starburst galaxies should produce very-high-energy (VHE;E > 100 GeV) gamma-ray emission via the interaction of their copious quantities of cosmic rays with the large reservoirs of dense gas within the galaxies. Such VHE emission was detected by VERITAS from the starburst galaxy M82 in 2008–09. An extensive, multiyear campaign followed these initial observations, yielding a total of 254 hr of good-quality VERITAS data on M82. Leveraging modern analysis techniques and the larger exposure, these VERITAS data show a more statistically significant VHE signal (∼6.5 standard deviations,σ). The corresponding photon spectrum is well fit by a power law (Γ = 2.3 ± 0.3_stat ± 0.2_sys), and the observed integral flux isF(>450 GeV) = (3.2 ± 0.6_stat ± 0.6_sys) × 10⁻¹³cm⁻²s⁻¹, or ∼0.4% of the Crab Nebula flux above the same energy threshold. The improved VERITAS measurements, when combined with various multiwavelength data, enable modeling of the underlying emission and transport processes. A purely leptonic scenario is found to be a poor representation of the gamma-ray spectral energy distribution (SED). A lepto-hadronic scenario with cosmic rays following a power-law spectrum in momentum (indexs ≃ 2.25) and with significant bremsstrahlung below 1 GeV provides a good match to the observed SED. The synchrotron emission from the secondary electrons indicates that efficient nonradiative losses of cosmic-ray electrons may be related to advective escape from the starburst core. 

Abstract The Accelerator Neutrino Neutron Interaction Experiment (ANNIE) is a 26-ton water Cherenkov neutrino detector installed on the Booster Neutrino Beam (BNB) at Fermilab. Its main physics goals are to perform a measurement of the neutron yield from neutrino-nucleus interactions, as well as a measurement of the charged-current cross section of muon neutrinos. An equally important focus is the research and development of new detector technologies and target media. Specifically, water-based liquid scintillator (WbLS) is of interest as a novel detector medium, as it allows for the simultaneous detection of Cherenkov light and scintillation. This paper presents the deployment of a 366 L WbLS vessel in ANNIE in March 2023 and the subsequent detection of both Cherenkov light and scintillation from the WbLS. This proof-of-concept allows for the future development of reconstruction and particle identification algorithms in ANNIE, as well as dedicated analyses within the WbLS volume, such as the search for neutral-current events and the hadronic scintillation component. 

Abstract G106.3+2.7, commonly considered to be a composite supernova remnant (SNR), is characterized by a boomerang-shaped pulsar wind nebula (PWN) and two distinct (“head” and “tail”) regions in the radio band. A discovery of very-high-energy gamma-ray emission (E_γ> 100 GeV) followed by the recent detection of ultrahigh-energy gamma-ray emission (E_γ> 100 TeV) from the tail region suggests that G106.3+2.7 is a PeVatron candidate. We present a comprehensive multiwavelength study of the Boomerang PWN (100″ around PSR J2229+6114) using archival radio and Chandra data obtained two decades ago, a new NuSTAR X-ray observation from 2020, and upper limits on gamma-ray fluxes obtained by Fermi-LAT and VERITAS observatories. The NuSTAR observation allowed us to detect a 51.67 ms spin period from the pulsar PSR J2229+6114 and the PWN emission characterized by a power-law model with Γ = 1.52 ± 0.06 up to 20 keV. Contrary to the previous radio study by Kothes et al., we prefer a much lower PWNB-field (B∼ 3μG) and larger distance (d∼ 8 kpc) based on (1) the nonvarying X-ray flux over the last two decades, (2) the energy-dependent X-ray size of the PWN resulting from synchrotron burn-off, and (3) the multiwavelength spectral energy distribution (SED) data. Our SED model suggests that the PWN is currently re-expanding after being compressed by the SNR reverse shock ∼1000 yr ago. In this case, the head region should be formed by GeV–TeV electrons injected earlier by the pulsar propagating into the low-density environment. 

Abstract In 2017 February, the blazar OJ 287 underwent a period of intense multiwavelength activity. It reached a new historic peak in the soft X-ray (0.3–10 keV) band, as measured by the Swift X-ray Telescope. This event coincides with a very-high-energy (VHE)γ-ray outburst that led VERITAS to detect emission above 100 GeV, with a detection significance of 10σ(from 2016 December 9 to 2017 March 31). The time-averaged VHEγ-ray spectrum was consistent with a soft power law (Γ = −3.81 ± 0.26) and an integral flux corresponding to ∼2.4% that of the Crab Nebula above the same energy. Contemporaneous data from multiple instruments across the electromagnetic spectrum reveal a complex flaring behavior, primarily in the soft X-ray and VHE bands. To investigate the possible origin of such an event, our study focuses on three distinct activity states: before, during, and after the 2017 February peak. The spectral energy distributions during these periods suggest the presence of at least two nonthermal emission zones, with the more compact one responsible for the observed flare. Broadband modeling results and observations of a new radio knot in the jet of OJ 287 in 2017 are consistent with a flare originating from a strong recollimation shock outside the radio core. 

Abstract Dark matter is a key piece of the current cosmological scenario, with weakly interacting massive particles (WIMPs) a leading dark matter candidate. WIMPs have not been detected in their conventional parameter space (100 GeV ≲ M χ ≲ 100 TeV), a mass range accessible with current Imaging Atmospheric Cherenkov Telescopes. As ultraheavy dark matter (UHDM; M χ ≳ 100 TeV) has been suggested as an underexplored alternative to the WIMP paradigm, we search for an indirect dark matter annihilation signal in a higher mass range (up to 30 PeV) with the VERITAS γ -ray observatory. With 216 hr of observations of four dwarf spheroidal galaxies, we perform an unbinned likelihood analysis. We find no evidence of a γ -ray signal from UHDM annihilation above the background fluctuation for any individual dwarf galaxy nor for a joint-fit analysis, and consequently constrain the velocity-weighted annihilation cross section of UHDM for dark matter particle masses between 1 TeV and 30 PeV. We additionally set constraints on the allowed radius of a composite UHDM particle. 

Abstract While the sources of the diffuse astrophysical neutrino flux detected by the IceCube Neutrino Observatory are still largely unknown, one of the promising methods to improve our understanding of them is investigating the potential temporal and spatial correlations between neutrino alerts and the electromagnetic radiation from blazars. We report on the multiwavelength target-of-opportunity observations of the blazar B3 2247+381, taken in response to an IceCube multiplet alert for a cluster of muon neutrino events compatible with the source location between 2022 May 20 and 2022 November 10. B3 2247+381 was not detected with VERITAS during this time period. The source was found to be in a low-flux state in the optical, ultraviolet, and gamma-ray bands for the time interval corresponding to the neutrino event, but was detected in the hard X-ray band with NuSTAR during this period. We find the multiwavelength spectral energy distribution is described well using a simple one-zone leptonic synchrotron self-Compton radiation model. Moreover, assuming the neutrinos originate from hadronic processes within the jet, the neutrino flux would be accompanied by a photon flux from the cascade emission, and the integrated photon flux required in such a case would significantly exceed the total multiwavelength fluxes and the VERITAS upper limits presented here. The lack of flaring activity observed with VERITAS, combined with the low multiwavelength flux levels, as well as the significance of the neutrino excess being at a 3σlevel (uncorrected for trials), makes B3 2247+381 an unlikely source of the IceCube multiplet. We conclude that the neutrino excess is likely a background fluctuation. 

Abstract The ground-based gamma-ray observatory Very Energetic Radiation Imaging Telescope Array System (VERITAS, https://veritas.sao.arizona.edu/ ) is sensitive to photons of astrophysical origin with energies in the range between ≈85 GeV and ≈30 TeV. The instrument consists of four 12 m diameter imaging Cherenkov telescopes operating at the Fred Lawrence Whipple Observatory in southern Arizona. VERITAS started four-telescope operations in 2007 and collects about 1100 hr of good-weather data per year. The VERITAS collaboration has published over 100 journal articles since 2008 reporting on gamma-ray observations of a large variety of objects: Galactic sources like supernova remnants, pulsar wind nebulae, and binary systems; extragalactic sources like star-forming galaxies, dwarf-spheroidal galaxies, and highly variable active galactic nuclei. This note presents VTSCat: the catalog of high-level data products from all VERITAS publications. 

Abstract We report the detection of very high energy gamma-ray emission from the blazar S3 1227+25 (VER J1230+253) with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). VERITAS observations of the source were triggered by the detection of a hard-spectrum GeV flare on 2015 May 15 with the Fermi-Large Area Telescope (LAT). A combined 5 hr VERITAS exposure on May 16 and 18 resulted in a strong 13σdetection with a differential photon spectral index, Γ = 3.8 ± 0.4, and a flux level at 9% of the Crab Nebula above 120 GeV. This also triggered target-of-opportunity observations with Swift, optical photometry, polarimetry, and radio measurements, also presented in this work, in addition to the VERITAS and Fermi-LAT data. A temporal analysis of the gamma-ray flux during this period finds evidence of a shortest variability timescale ofτ_obs= 6.2 ± 0.9 hr, indicating emission from compact regions within the jet, and the combined gamma-ray spectrum shows no strong evidence of a spectral cutoff. An investigation into correlations between the multiwavelength observations found evidence of optical and gamma-ray correlations, suggesting a single-zone model of emission. Finally, the multiwavelength spectral energy distribution is well described by a simple one-zone leptonic synchrotron self-Compton radiation model.